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Abstract 

Let there be given a probability measure on the unit circle T 
of the complex plane and consider the inner product induced by /U. 
In this paper we consider the problem of orthogonalizing a sequence 
of monomials {z'^''}k, for a certain order of the € Z, by means of 
the Gram-Schmidt orthogonalization process. This leads to a basis 
of orthonormal Laurent polynomials {^Jk}k- Wc show that the ma- 
trix representation with respect to the basis {tpk}k of the operator of 
multiplication by z is an infinite unitary or isometric matrix allowing 
a 'snake-shaped' matrix factorization. Here the 'snake shape' of the 
factorization is to be understood in terms of its graphical representa- 
tion via sequences of little line segments, following an earlier work of 
S. Delvaux and M. Van Barel. Wc show that the shape of the snake is 
determined by the order in which the monomials {z'^''}k are orthogo- 
nalized, while the 'segments' of the snake are canonically determined 
in terms of the Schur parameters for /i. Isometric Hesscnberg matrices 
and unitary five-diagonal matrices (CMV matrices) follow as a special 
case of the presented formalism. 
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1 Introduction 



1.1 Unitary Hessenberg and five-diagonal matrices 

In recent years, there has been a lot of research activity on the topic of 
unitary five- diagonal matrices, also known as CM V matrices. These matrices 
have been used by researchers in various contexts, see e.g. [6], [Z]-l8], [16] . 
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Explicitly, the CMV matrix looks like 
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where a^, fc = 0, 1,2, . . . are complex numbers satisfying \ak\ < 1 (the so- 
called Schur parameters or Verblunsky coefficients) and pk ■= ^/Y^-Jc^k\^ S 
(0, 1] are the so-called complementary Schur parameters. The matrix C = 
(cij)jj>c[^ in ([1]) can be seen to be unitary and five-diagonal, in the sense 
that Cij = whenever \i — j\ > 2. More precisely, the nonzero entries of C 
follow a kind of zigzag shape around the main diagonal. 

The terminology 'CMV matrix' for the matrix in ([T]) originates from the 
book of Simon [27], who named these matrices after a 2003 paper by Can- 
tero. Moral and Velazquez [7]. But this terminology is far from historically 
correct, since the latter paper [7] is in fact a rediscovery of facts which were 
already known by the numerical analysis community in the early 1990's; a 
survey of these early results can be found in the review paper by Watkins 
|33] : see also p8]. 

In the present paper, we prefer to avoid such historical discussions and 
we will therefore use the neutral term 'unitary five-diagonal matrix' to refer 
to these CMV matrices. 

Unitary five-diagonal matrices have a number of interesting features, in- 
cluding the statement proven in the literature that (see further in this paper 
for more details) from all non-trivial classes of unitary matrices, unitary five- 
diagonals have the smallest bandwidth. Here the word 'non-trivial' refers to 
matrices which are not expressible as a direct sum of smaller matrices. 

While this statement about the minimal bandwidth is certainly correct, 

*In the rest of the paper and for convenience with the notation, we will label the rows 
and columns of any matrix starting with index 0. As an example, the element ci,i in the 
matrix will take the value — aool. 
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it is a curious fact that this does not imply that unitary five-diagonal ma- 
trices are also numerically superior with respect to other non-trivial classes 
of unitary matrices. For example, another class of unitary matrices which is 
often used in the literature is the class of unitary Hessenberg matrices, given 
explicitly by 



Note that the matrix in Q is of infinite dimension. Therefore it should 
actually be called an isometric Hessenberg matrix rather than a unitary 
Hessenberg matrix, see e.g. |18j . Nevertheless, throughout this paper we 
choose to slightly abuse the notation by using everywhere the unified term 
unitary rather than isometric, for reasons of notational uniformity. 

In (§ we use again the notation q^, to denote the Schur parameters 
and complementary Schur parameters, respectively. These are the same 
numbers as in the matrix ([T|; see further. 



whenever i — j > 2. Note however that the upper triangular part of this 
matrix is in general dense. 

Now the point is that unitary Hessenberg matrices as in (j2]) are known 
to be just as efficient to manipulate as unitary five-diagonal matrices! Al- 
though this fact is known by numerical specialists, it seems that it is not so 
well-known in part of the theoretical community. Therefore, let us describe 
this now in somewhat more detail. 

The naive idea would be that unitary Hessenberg matrices are 'inefficient' 
to work with since these matrices have a 'full' upper triangular part, in 
contrast to unitary five-diagonal matrices. But this would be a too quick 
conclusion. Having a better look at the problem, one can note that the 
upper triangular part of a unitary Hessenberg matrix is rank structured in 
the sense that each submatrix that can be taken out of the upper triangular 
part of such a matrix, has rank at most equal to 1. This can be easily 
verified using e.g. the explicit expressions of the entries of the matrix Ti. in 



Going one step further, one can note that the rank structure in the upper 
triangular part of TC is in fact a consequence of an even more structural 
theorem. Denote with Gk,k+i ^ Givens transformation (also called Jacobi 
transformation) 
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where the / denote identity matrices of suitable sizes (the second one being 
infinite), and where Gk,k+i is a 2x2 unitary matrix positioned in rows and 
columns {k, k+ 1}. Thus the matrix Gk^k+i differs from the identity matrix 
only by its entries in rows and columns {fc, /c + 1}. Givens transformations 
can be considered as the most elementary type of unitary matrices. They 
can be used as building blocks to construct more general unitary matrices. 
Of interest for the present discussion is the fact that any (infinite) unitary 
Hessenberg matrix TL allows a factorization as a product of Givens transfor- 
mations in the form 

W = Go,lGl,2G2,3G'3,4 .... (4) 

This result can be shown using only some basic linear algebra |15| [T7] . 

Applying this factorization to the matrix in (|2|, one can actually 
specify this result by noting that the feth Givens transformation Gk,k+i in 
Q must have nontrivial part given by 

(5) 

In other words, the 'cosines' and 'sines' of the Givens transformations in (Q 
are nothing but the Schur parameters and complementary Schur parameters, 
respectively. This result was first established in the present context by 
Ammar, Gragg and Reichel 

Incidently, note that the Givens transformations in ^ are of a special 
form in the sense that they have real positive off-diagonal elements and 
determinant —1. 

We also note the following finite dimensional equivalent of (|4]): any uni- 
tary Hessenberg matrix Ti of size n x n allows a factorization in the form 

W = G'o,lGi,2G2,3G'3,4 • • • Gn-2,n-l, (6) 

for suitable Givens transformations G^^k+i: k = 0,1, . . . ,n — 2. 

The main point of ^ is that it shows that unitary Hessenberg matrices 
of size n can be compactly represented using only 0{n) parameters, just as 
is the case for unitary five-diagonal ones. Working with such an 0(n) ma- 
trix representation, the eigenvalue problem for unitary Hessenberg matrices 
can be solved numerically in a fast and accurate way; see the end of Section 
|4] for some references to eigenvalue computation algorithms in the litera- 
ture. These algorithms can be canonically expressed in terms of the matrix 
factorization (|6j), i.e., in terms of the Schur parameters of the problem. 



G 



k,k+l 



ak 
Pk 



Pk 
-Q-k 



1.2 Graphical representation 

In [13] , a graphical notation was introduced where matrix factorizations with 
Givens transformations are represented via sequences of little line segments. 
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The graphical representation is obtained as follows. Let A be some 
arbitrary matrix (which will play no role in what follows) and suppose that 
we update A i— > Gk,k+iA. This means that the A:th and {k + l)th row of A 
are replaced by linear combinations thereof, while the other rows of A are 
left unaltered. We can visualize this operation by drawing a vertical line 
segment on the left of the two modified rows of A. 

One can then apply this idea in an iterative way. For example, when 
updating A by means of an operation A i-^ Gk+i,k+2Gk,k+iA, one places first 
a vertical line segment on the left of rows k,k + l (this deals with the update 
A i-H- Gk,k+iA), and subsequently places a second vertical line segment on 
the left of the former one, this time at the height of rows k + l,k + 2. We 
obtain in this way two successive vertical line segments. Clearly, any number 
of Givens transformations can be represented in such a way. 

Now the key point is that we identify each Gk^k+i with its corresponding 
vertical line segment. We hereby make abstraction of the matrix A on whose 
rows these operations were assumed to act. For example, the graphical 
representation of the factorization ^ with n = 8 is shown in Figure [T| 



1 I 

2 I 

3 I 

4 I 

5 I 

6 I 
7 



Figure 1: The figure shows in a graphical way the decomposition as a product 
of Givens transformations of the unitary Hessenberg matrix in ([g]) with 
n = 8. 

Concerning Figure [T] note that the top leftmost line segment in this 
figure (which is assumed to be placed at 'height' and 1; cf. the indices 
on the left of the figure) corresponds to the leftmost factor Go,i in ([s]). 
Similarly, the second line segment corresponds to the factor Gi,2 in ([6|, and 
so on. We again emphasize that the line segments in Figure [T] should be 
imagined as 'acting' on the rows of some (invisible) matrix A. See |13l [H] 
for more applications of this graphical notation. 

It is known that also unitary five-diagonal matrices allow a factorization 
as a product of Givens transformations. More precisely [U HI [3 [33], the 
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matrix C in ([T| allows the factorization 
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which can be rewritten as 

C = {. . . ^6,76*4,502,3^0,1) • (Gi, 2^3,465,6 . . .)) (7) 

where the Gk,k+i are again defined by ([s]) and ([5]). The factorization ([7| is 
represented graphically for n = 8 in Figure [2j 



(a) 



Figure 2: The figure shows in a graphical way (a) the decomposition as a 
product of Givens transformations of the unitary five-diagonal matrix ([T]), 
(b) the 'snake shape' underlying this decomposition. 

Let us comment on Figure [2j The leftmost series of line segments in 
Figure 2(a) corresponds to the leftmost factor in ([7]). The order in which 



these Givens transformations are multiplied is clearly irrelevant; therefore we 
are allowed to place them all graphically aligned along the same vertical line. 
Similarly, the rightmost series of line segments in Figure |2(a)| corresponds to 
the rightmost factor in ([7| . To explain Figure 2(b)[ imagine that one moves 



from the top to the bottom of the graphical representation. Then one can 
imagine a certain zigzag 'snake shape' underlying the factorization, which 
is shown in Figure 2(b)| 



Note that in the above discussions we did not describe the way how uni- 
tary Hessenberg and unitary five-diagonal matrices arise in practice as ma- 
trix representations of a certain operator. At present, it will suffice to know 
that they are matrix representations of the operator of multiplication by z, 
expressed in a certain basis of orthonormal Laurent polynomials {ipkiz)}k- 
This orthonormal basi^is obtained by applying the Gram-Schmidt orthog- 
onalization process to the sequence of monomials 

1, z, z'^ , , . . . , and 1, z, z^'^ , z'^ , z^'^ , . . . , (8) 



^When using the word 'basis' it should of course be said of which vector space this 
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for the unitary Hessenberg and unitary five-diagonal case, respectively. 



1.3 Snake-shaped matrix factorizations 

The aim of this paper it to carry the above observations one step further. 
We will show that with respect to a general orthonormal basis of Laurent 
polynomials {ipkiz)}ki obtained by orthogonalizing a general sequence of 
monomials (satisfying some conditions to be described in detail in Section 



2.1 ), the operator of multiplication by z becomes an infinite unitary matrix 
allowing a snake-shaped matrix factorization. We will use the latter term 
to denote an infinite matrix product S = IlfcLo^fc.fc+i' where the factors 
under the f^-symbol are multiplied in a certain order. Here the 'segments' 
Gk^k+i of the snake are canonically fixed in terms of the Schur parameters 
by means of ([3| and while the 'shape' of the snake, i.e., the order in 
which the Gk,k+i are multiplied, will be determined by the order in which 
the monomials have been orthogonalized. 

To fix the ideas, consider the sequence of monomials 

1, z~\ z, z^, z\ z-\ z-\ z\ z^ . . . . (9) 

With respect to the resulting orthonormal basis of Laurent polynomials 



{'4>k{z)}k (see Section 2.1 for details), the operator of multiplication by z 
will then allow a snake-shaped matrix factorization S = S^°°\ We claim 
that this factorization is built by means of the following recipe: 

1. Considering the monomial 1 = z*^ in the position of ([9|, we initialize 
S^^^ := Go,i. Then we apply the following procedure for k > 1: 

2. If the /cth monomial in ([9| has a positive exponent, we multiply the 
matrix with a new Givens transformation on the right by setting 
Sik) :=s(^-^)Gk,k+i; 

3. If the /cth monomial in ^ has a negative exponent, we multiply the 
matrix with a new Givens transformation on the left by setting S^''^ := 
Gk,k+iS^''~^^ ■ 

For the sequence of monomials Q , this recipe gives rise to the following 
series of iterate matrices S^'^^: 

S^^^ = Go,i, 5*^^^ = Gi^2 ■ Co,!, 

5*^^) = Gi^2 ■ Go, 1^2,3, S^^^ = G3,4Gi,2 • Go,lG2,3 

5(4) = ^3,4^1,2 • Go,lG2,3G4,5! S^^^ = ^3,4^1,2 • Go,lG2,3G4,5G5,6, 

5(6) = 06,7^3,4^1,2 • Go,lG2,3G4,5G5,6, • • . . 

basis is. The obvious vector space here is the space of quadratically integrable functions 
with respect to the probabihty measure fi to be considered in Section [2. 1[ However, this 
leads to problems since e.g. the sequence of monomials 1, 2, 0^, ... is in general not a basis 
for this vector space! (see e.g. [30]). Nevertheless, in this paper we will not be concerned 
about such topics, and we will therefore use the word 'basis' in a somewhat extended way. 
We hope that this will not lead to inconvenience for the more scrutineer reader. 
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This leads to the final matrix factorization 

S = S^°°^ = (• • • ^7,8^6,703,401^2) • (00,102,304,505,608,909,10 • • •)• (10) 
This factorization is shown graphically in Figure [3| 




(a) 



Figure 3: The figure shows in a graphical way (a) the decomposition as a 



product of Givens transformations of the matrix S in (10), (b) the 'snake 
shape' underlying this decomposition. 

Let us comment on Figure [3} The 'snake' in this figure was built by 
means of the following recipe: 

1. Starting with a snake consisting of a single line segment Oo,i, we apply 
the following procedure for A: > 1: 

2. If the A;th monomial in (joj) has a positive exponent, the snake moves 
towards the bottom right, i.e., we add a new line segment on the bottom 
right of the snake; 

3. If the kth monomial in ([9| has a negative exponent, the snake moves 
towards the bottom left, i.e., we add a new line segment on the bottom 
left of the snake. 

Of course, this recipe is nothing but a direct translation of the recipe that 



led us to the matrix factorization (10). 



The reader should check that the above procedures are also valid for the 
unitary Hessenberg and the unitary five-diagonal case (cf. (|8| and Figures 
[Tl[2l). 



1.4 Outline of the paper 



The fact that the recipe in Section 1.3 leads to the correct matrix repre- 



sentation of the operator of multiplication by z with respect to the basis 
of orthonormal Laurent polynomials {'4'k{z)}k will be shown in Section [2j 
Our proof makes use of essentially three facts: (i) an observation of Cruz- 
Barroso et al. [12, (see also Watkins ^33j) expressing the intimate connection 
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between orthonormal Laurent polynomials and Szego polynomials; (ii) the 
well-known Szego recursion [30]; and (iii) an argument of Simon |28] us- 
ing 'intermediary bases' in the unitary Hessenberg case. The full proof is 
however rather technical and requires some administrational book-keeping. 

By factoring out a snake-shaped matrix product like ( 10 ), one can obtain 
explicit expressions for the entries of the matrix, generalizing the expansions 
in ([T| and ^ . This will be the topic of Section [s] where we will describe 
a graphical rule for determining the zero pattern of the matrix S as well as 
the shape of its non-zero elements. 

Finally, in Section |4] we will briefly consider some connections between 
snake-shaped matrix factorizations and Szego quadrature formulas. We will 
show that the known results involving unitary Hessenberg and unitary five- 
diagonal matrices can all be formulated in terms of a general snake-shaped 
matrix factorization S, extending an observation of Ammar, Gragg and 
Reichel [I]. 

The remainder of this paper is organized as follows. Section [2] discusses 
some preliminaries about sequences of orthogonal Laurent polynomials on 
the unit circle and proves the main result about snake-shaped matrix fac- 
torizations. Section [3] discusses the entry-wise expansion of snake-shaped 
matrix factorizations. Finally, Section |4] considers the connection with Szego 
quadrature formulas. 



2 Snake-shaped matrix factorizations: main result 

This section is devoted to the proof of our main result about snake-shaped 
matrix factorizations, showing how these occur as the matrix representation 
of the operator of multiplication by z with respect to a basis of orthonormal 
Laurent polynomials. We start with some preliminaries. 



2.1 Sequences of orthogonal Laurent polynomials on the unit 
circle 

In this first subsection we fix some notations and conventions concerning 
orthogonal Laurent polynomials on the unit circle (see [6], jl0j-|12]). We 
denote by T := {z G C : \z\ = 1} the unit circle in the complex plane 
and by A := C[z,z^^] the complex vector space of Laurent polynomials 
in the variable z. For a given order n £ N and an ordinary polynomial 
P{^) — Ylk=o'^kz'' , we define its dual as p*{z) := z"'p{l/z), or explicitly 
P*{z) = ^^=0 '^n-kz^- Here the bar denotes complex conjugation. 

Throughout the paper, we shall be dealing with a finite positive non- 
discrete Borel measure ^ supported on the unit circle T (which induces a 
measure on the interval [— vr, vr] that we also denote by /i), normalized by 
the condition J^^ d^[6) = 1, i.e., a probability measure. As usual, the inner 
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product induced by fi is given by 



if, 9)= / f{e^')g{e^')dfi{9). (11) 



For our purposes, we start constructing a sequence of subspaces of Lau- 
rent polynomials {£„}^o satisfying 

Cq := span{l} , dim (£„) = n + 1 , C Cn+i , n > 1. 

This can be done, by taking a sequence {pn}^=o of non-negative integers 
such that Po = 0, < p„ < n and s„ := pn — Pn-i £ {0, 1} for all n > 1. In 
the sequel, a sequence {pn}^=o satisfying these requirements will be called a 
generating sequence. Observe that in this case both {pn}^=o and {n—pn}'^=Q 
are non-negative non-decreasing sequences. Then, set 

Cn := span {z^ : -p„ <j<n- _p„} 

and set £_i := {0} to be the trivial subspace. Observe that A = U^o 
if and only if lim pn = lim {n — pn) = oo and that for all n > 1, 



r 



Cn-i e span{z"-P"} if s„ = 0, 
£„_i © span{z"P"} if s„ = 1. 



By applying the Gram-Schmidt orthogonalization procedure to an 
orthonormal basis {%1)q{z), . . . ,ipn{z)} can be obtained. If we repeat the 
process for each n > 0, a sequence {V'n(-z)}^o Laurent polynomials can 
be obtained satisfying, for all n, m > 0: 



1. i>n{z) G Cn \ Cn~l, tpo{z) = 1, 

2. ipniz) has a real positive coefficient for the power 



if s„ = 
z-P" if s„ = 1 ' 



3. {'4>n{z),^pm{z)) 



if n 7^ m 

1 \i n = m 



This sequence will be called a sequence of orthonormal Laurent polyno- 
mials for the measure n and the generating sequence {pn}^=o- 
Let us illustrate these ideas with three examples. 

Example 1 Consider the sequence of monomials given by ^ and the mono- 
mial 1 = z^ in the position 0. Then, the construction of the sequence {sn}'^=i 
is nothing but to take Sn = if the nth monomial in ^ has a positive expo- 
nent and s„ = 1 if it is negative, whereas pn counts the number of negative 
monomials positioned up to n. Hence, {sn}'^=i = {1, 0, 1, 0, 0, 1, 1,0, 0, . . .} 
and {pn}'^=o = {0, 1, 1,2, 2, 2, 3, 4, 4, 4, . . .}. 
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Example 2 If = for all k > 1, then Cn is the space of ordinary 
polynomials of degree at most n. In this case the Gram-Schmidt orthogonal- 
ization process is applied to the sequence of monomials {1, z, z'^, z^, . . .} and 
the resulting orthonormal Laurent polynomials ipn{z) are just the well-known 
orthonormal Szego polynomials ipn{z); see e.g. ISO]. 

Example 3 // Sfc = A; + 1 mod 2 for all k > 1, then the Gram-Schmidt 
orthogonalization process is applied to the sequence {1, z, z~^ , z"^ , z~'^ , . . .} , 
where the monomials z^ and z^^ occur in an alternating way. The resulting 
sequence {tpn{z)}^^Q was firstly considered by Thron in f3lf and it is called 
the CMV basis in The CMV basis can actually be expressed in terms 

of the Szego polynomials as (see e.g. Q El US Ell IM/) 



In the general case, one has the following result. 

Lemma 4 (Cruz-Barroso et al. Ill 2^ : see also Watkins ]33^) The family 
{'0n(-z)}5^o sequence of orthonormal Laurent polynomials on the unit 

circle for a measure fi and the ordering induced by the generating sequence 
{Pn}n=0' '^^d only if 



{(/9„(z)}^Q being the sequence of orthonormal Szego polynomials for fi. 



Lemma|4]shows that the orthonormal Laurent polynomials {ipn{z)}n are 
very closely related to the usual Szego polynomials {(/?n(-z)}n and their duals, 
and this for any choice of the generating sequence {pnj^o- ^^^^ need 
this result in what follows. 

2.2 The main result 

In this subsection we state and prove the main result of this paper. Let 
{ipn{z)}'^^Q be the sequence of orthonormal Laurent polynomials on the unit 
circle for the measure /i and the ordering induced by the generating sequence 
{Pn}'^=Q- To distinguish them from the other bases to be constructed in 
this section, we will equip these Laurent polynomials with a superscript: 
'ipn\z) := Tpniz). We will also find it convenient to use the vectorial notation 
'ip^^\z) := {'4^n\z))'^^Q. Thus, tp^^^ is an infinite dimensional vector whose 
nth component is the nth orthonormal Laurent polynomial ^pn^ (n > 0). 



ipo{z),ipi{z),z ^ip*2{z),z ^ips{z),z ^ipl{z),z ^V?5(2;),.... 




(12) 



□ 
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Let M denote the operator of multiplication by z on the space of quadrat- 



ically integrable functions with respect to the inner product (11). Thus M 
is defined by the action 

M : f{z) ^ zf{z) , /GL^(T). 

Since we are working on the unit circle T, the operator M is actually unitary. 

Let S be the matrix representation of the operator M with respect to 
the orthonormal basis i.e., 



[{4'\z),z^f\z))] 



oo 

i,j=0 



(13) 



where the inner product is defined in ( |11[ ). Here the expression on the 
second line should be regarded as a compact vectorial notation of the line 
above. Note that 5 is a unitarjj^] matrix: this follows since it is the matrix 
representation of the unitary operator M with respect to the orthonormal 
basis tp^^\ 

Now we are in position to prove the following result. We will do this by 
using a modification of an argument of Simon [28| for the unitary Hessen- 
berg case. The main ingredient of the proof will be the well-known Szego 
recursion, expressed in the form (see e.g. j30j ) 



Z(pk{z) 






Pk 










Pk 


-ak _ 




_ Vk+iiz) 



(14) 



where <fk{z) and ^l.{z) denote the orthonormal Szego polynomial of degree k 
and its dual respectively. Note that the coefficient matrix in ( 14 ) is nothing 
but the nontrivial part ([5| of the Givens transformation Gk^k+i- 

Theorem 5 Let {^n(-2)}^o sequence of orthonormal Laurent poly- 

nomials on the unit circle for a measure n and the ordering induced by 
the generating sequence {pn}^=Q- Then the matrix S in (13) can be factored 
into a snake-shaped matrix factorization S = S^°°\ constructed by the recipe 
given in Section \l.3\ 

PROOF. We will construct a sequence of intermediary bases tl^^''\ k > 1, 
in such a way that for each k, there exists an index I S {0, 1, . . . , A; — 1} such 
that tf^^^^ is the same as except for a change in the {k — l)th and kth 

*We should be careful here about the word 'unitary'. If A 7^ UJJ^o then ipl^^ will in 
general not be a dense subspace of (IT)- In that case the column space of the matrix S 
does not correspond to the full space (IT), and hence the matrix 5 should then be called 
isometric rather than unitary, as e.g. in [TH]. See also the paragraph following Eq. ([2|. 
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components. These intermediary bases will serve to factorize the matrix S. 
For example, note that (13) can be rewritten as 



{iP^^\z),iP^^\z)) ■ {ip^^\z),zip'^^\z)), 



(15) 
(16) 



for any choice of the basis Indeed, the jth column of the matrix (16) 



is obtained by expressing the function zil^j (z) in terms of the basis 
which is then in its turn expressed in terms of the basis tl)^^\z). Obviously 
this gives the same result as directly expressing zipj^\z) in terms of the basis 

%p^^\z), i.e., it equals the jth column of (15) (we recall again our convention 
with the notation: j > 0). 



Note that instead of (16) we could also have written a slightly modified 
version of it: 



(^/>(0)(z),z^(0)(z)) 

(^/)(°)(z),z^W(z)) • {zip'^^\z),ztp'-°\z)) 
(^(°)(z),z^W(z)) • (V^W(z),^(°)(z)) 



(17) 



for any choice of the basis i/^*-^^ and where we have used the general fact that 
{zf{z),zg{z)) = {f{z),g{z)) for any functions /, 5 : T — > C, which follows 



from (11) and the fact that z G T. The choice between (16) and (17) will 
depend on the fact whether si = or si = 1, respectively; see further. 

The point will now be to make a good choice for the intermediary bases 
For example, the 'good' choice for •0^^-' will be the one for which one 



of the factors in (16) (or (17)) equals the Givens transformation Go,i, while 



the other factor is of the form 



1 

* 



where * denotes an irrelevant submatrix (which is actually of infinite di- 



mension). Explicitly, the basis ip 



(1) 



is given by tp^^ 



A-2s 



(1) 



ijj^^'^ for all A; > 2. Repeating this idea inductively 
for all subsequent bases i/'^'^^ will ultimately lead to the decomposition of S 
as an infinite product of Givens transformations. 

Let us now formalize these ideas. We work with the induction hypothesis 
that after the /cth step, /c > CR we have decomposed the matrix S as 



5 = 5(^"i)xW, ifsfc 
5 = xW5(^-i), ifsfc 



0, 
1, 



^This procedure also works for fc = provided that we set 5' 
So = 1, since either choice will give the same result. 



(18) 



/ and So = or 
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where iS^'^^^) is the (A; — l)th iterate matrix of t he s nake-shaped matrix 
factorization S^°°^ (cf. the construction in Section 1.3), while X^'^^ equals 
the identity matrix in its first k x k block, i.e.. 



Ik 
* 



(19) 



with /fc the identity matrix of size k. We also assume by induction that 

X^^^ = {^p^^\z),z^^"'\z)), (20) 

where I, m are certain indices in {0, 1, . . . , A;} with at least one of them equal 
to k (we could actually give explicit expressions for /,m but will not need 



these in what follows) . Note that by combining the hypotheses ( 19 ) and ( 20 ) , 



we deduce that ipf'^ = zij:)^'^' {z) for alH G {0, 1, . . . , A; — 1}. In addition, we 
have the following induction hypothesis on the kth. components of ^/''•'^ and 



, (m) 



V'f(^) 



(21) 



where ipk{z) denotes the orthonormal Szego polynomial of degree k. 

Given all the above induction hypotheses, we can now come to the in- 
duction step A; 1-^ A; -|- 1. To this end, we should try to peel off a new Givens 
transformation Gk,k+i from the matrix 5. Assume that the first k interme- 
diary bases ■ ■ ■ , have already been constructed. We want to define 
the next intermediary basis i/j^'^'*'^^. We distinguish between two cases: 

1. If Sfc+i = 0, we define to be the same as ip^^\ except for its A:th 

and (A; -|- l)th components. More precisely, we set 



(3) 



G 



k,k+l 



V'f(^) 



Gk,k+i • z-P>'+' 

2;-Pk + l 



Zipk{z) 



(22) 

(23) 
(24) 



Here the second equality follows from the first lines of (21) and (12) 



(recall that the (A: -|- l)th component of has not been changed yet 
with respect to 'ip^^^), and from the fact that Pk+i = Pk by assump- 
tion. On the other hand, the third equality is nothing but the Szego 



recursion (14). 
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Then, we can factorize (20) as 



(25) 



where the third equahty follows from ( 22 ) , and where the matrix 



^C^+i) In iiiQ fourth equality now equals the identity matrix in its 
first (A; + 1) x (A; + 1) block. The latter follows by the induction hy- 
pothesis in the first k x k block (rows and columns to A: — 1) and 



from the fact that, by the first line of (24) and the second line of (pT 
we have 



implying that also the feth column of the matrix X^^^^^ has all its 
entries equal to zero, except for the diagonal entry which equals one. 
By the unitarity, also the A;th row has then all its entries equal to zero, 
except for the diagonal entry. 

We can then replace the index / by its new value A;-|-l. We have already 



checked that the induction hypotheses (19) and (20) are inherited in 
this way as A; i-^ A; + 1 . Also the hypothesis (21 ) can be easily checked 



to remain valid in this way, by virtue of the second line of ( 24 ) and 



the first line of (12). Finally, we have to check that (18) remains also 
valid. To prove this, we use (18), (25), the construction of S^^^ in 



Section 1.3 and we distinguish between two cases: 



(a) If Sfc = then 



^(fc-i)j^(fc) 



(b) If Sfc = 1 then 



5 = ^('^^^(^-i) 
=: 

where we have used the commutativity of 5'^^^"'^) and X^'^^^^ since 
these matrices have a complementary zero pattern (the former 
equals the identity matrix except for its first (A;-|-l) x (A;-|-l) 
block, while the latter is precisely the identity matrix there, cf. 
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2. If Sfc+i = 1, we define ^(''+^) to be the same as except for its 

kth and {k + l)tli components. More precisely, we set 



G 



-1 

k,k+l 



Z(/7fc(z) 



(26) 

(27) 
(28) 



where we have used the second hnes of (21) and (12) (recall that the 
{k + l)th component of ■0^'"^ has not been changed yet with respect 
to tp^^^), the fact that Pk+i = Pfc + 1 by assumption and the Szego 



recursion (14). 



Then, we can factorize (20) as 



(^«(z),z^(™)(^)) 
(-0(')(z),z7/j('^+i)(z)) • {ziP^''+^\z),zip^"'\z)) 
{iP^^\z),zijj'^''+^\z)) • (-0('=+i)(z),-(/'^"H-z)) 
(^W(z),z^('=+i)(z)).Gfe,fc+i 



where the fourth step follows from (26). 



It is easy to check again that X^^"*"^^ equals the identity matrix in its 
first {k + 1) X (A; + 1) block by using the induction hypothesis in the 



first k X k block and from the first lines of (28) and (21) for the A:th 
row and column. 

We can then replace the index m by its new value k + 1. If follows 



from the above discussion that the induction hypotheses (19) and (20) 
are inherited in this way as k k + 1. Also the hypothesis (21 ) goes 



through, by virtue of the second lines of (12) and (28). Finally, the 



proof that also ( 18 ) goes through can be proven in a completely similar 



way as in the previous case. 

We have now completely established the induction hypothesis k 
hereby ending the proof of Theorem |5] 



k + 1, 
□ 



3 Entry-wise expansion of a snake-shaped matrix 
factorization 

In this section we discuss the entry- wise expansion of a snake-shaped matrix 
factorization S, hereby generalizing the expansions in ([T]) and ([2]). 



16 



3.1 Graphical rule for the entry-wise expansion of S 



First we will present a graphical rule for predicting both the position and 
the form of the non-zero entries of a snake-shaped matrix factorization S. 



We will illustrate the ideas for the matrix S given by (10) and Figure |3j 
A straightforward computation shows that the full expansion of this matrix 
S is given by (compare with [6], Example 4.5) 





"0 


Po 




















Pool 


—aoOi 


Pia2 


P1P2 
















PoPi 


-aoPi 


—aio^ 


-0LXP2 






















P201Z 


—a-ioz 




PiPiO^ 


p-iPiPb 













P2P3 


-Oi2PZ 




-asPioE 


-aaPiPb 



















Pi 




-a4P5 






















PbO^ 




P6 



















P^Peor 























P5P6P7 




-aep7 



V 



\ 



(29) 

Now the attentive reader will notice that the zero pattern of this matrix 
S has some similarity with the shape of its underlying snake as shown in 
Figure [3] Actually, we claim that the entry of the matrix S can be 
obtained from the following recipe (the 'E' stands for 'entry-wise'): 

El. Draw the snake underlying the matrix S (cf. Figure [3|; 

E2. Place a right-pointing arrow on the left of the snake at height i; 

E3. Place a left-pointing arrow on the right of the snake at height j; 

E4. Draw the path on the snake induced between these two arrows; 

E5. If the path moves monotonically from left to right, then the (i, j) entry 
of S equals a product of entries of the encountered Givens transfor- 
mations 

'" Pk 

Pk -ak 

on the path (see Step E5' below for a specification of this rule); 



Gk,k 



+1 



(30) 



E6. If the path does not move monotonically from left to right, then the 
entry of S equals zero. 

Let us illustrate this recipe for the (7, 5) entry of the matrix S (recall 
that we label the rows and columns of this matrix starting from the index 
0). The recipe is shown for this case in Figure |4j 

Let us comment on Figure E Figure 4(a) shows the snake shape of the 
matrix S (compare with Figure 3|, corresponding to Step El in the above 
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1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



(a) 




(b) 



(c) 



Figure 4: The figure shows (a) the snake shape underlying the matrix S in 



Equation (10), (b) the arrows on the left and on the right of the snake at 



height 7 and 5, respectively, and (c) the path on the snake induced between 
these two arrows. From this information, the value of the (7, 5) entry of the 
matrix S can be determined. 



recipe. Figure |4(b)| shows the arrows on the left and on the right of the 
snake at height 7 and 5, respectively, corresponding to Steps E2 and E3. 
The path on the snake induced between these two arrows is shown in Figure 
c) , corresponding to Step E4. Note that this path moves monotonically 
from left to right and passes through the Givens transformations 6*7^81 Gqj 
and G^fi. From Step E5 it then follows that the (7,5) entry of the matrix 
5 is a product of entries of these three Givens transformations. Actually, it 



equals p^p^aj (compare with (29)). 

As a second example, let us consider the (7,4) entry of the matrix S. 
The recipe is illustrated for this case in Figure [5] 




(a) 



(b) 



Figure 5: For the matrix S in Equation (|3]), the figure shows (a) the arrows 
on the left and on the right of the snake at height 7 and 4, respectively, and 
(b) the path on the snake induced between these two arrows. Since now 
the path does not move monotonically from left to right, it follows that the 
(7,4) entry of S equals zero. This corresponds again with (29). 
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Note that in the above example concerning the (7, 5) entry of the matrix 
S, we noticed that this entry equals the product of the (complex conjugate 
of) the Schur parameter cSy, on the one hand, and the complementary Schur 
parameters pQ, p5, on the other hand. To complete our description, let us 



now state an a priori rule to determine which of the four entries in (30) each 
Givens transformation Gk,k+i on the path in Step E5 contributes. 

Let us explain this rule for the first Givens transformation ^7^3 through 



which the path in Figure 4(c) passes (note that Gj^g corresponds to the 



bottom leftmost line segment on the path in Figure 4(c)). First, we will 
determine the row index of the entry contributed by To this end, 

imagine that we are in the line segment corresponding to G7.8 and that we 



move leftwards on the path. It is then seen from Figure 4(c) that we leave 
this line segment through its topmost index; hence we claim that the sought 
entry of ^7^8 will be in its topmost row. 

Next, to find the column index of the entry contributed by G73, imagine 
again that we start in the line segment corresponding to ^7^8 but move this 



time rightwards on the path. Since the path in Figure [4(c) proceeds upwards 



from left to right, we move then to the position of smaller indices. Hence 
the sought entry of G7,8 will be in its column with the smallest index, which 
is column 0. We conclude that the sought entry of ^7^3 lies in the (0, 0) 



position of ( 30 ) ; this gives us 07 . 

The entries contributed by Gqj and G^fi can be found in a similar way. 
The reader can check that in both cases, the relevant entries of Ggj, ^5,6 



are positioned in the (1,0) entry of (30) 



To summarize these ideas, let us introduce some notations. Denote with 
Gr,r+i and Gt,t+i the two outermost line segments of the path in Step E5. 
Note that r £ {i — l,i} and t G {j — 1, j}, with the precise value of r and 
t depending on the shape of the snake. For the example of Figure 4(c) we 
have r = i = 7 and t = j = 5. 

Denote with IC the set of indices k of the innermost Givens transforma- 
tions Gk,k+i on the path. Explicitly, /C equals {r-|-l,...,t — 1} ifr < t and 
+ 1, . . . , r — 1} if r > t (it is understood that /C = when |r — t| =1). 

It is easily seen from the above graphical rule that the {Gk,k+i}keK. in 
Step E5 always contribute their complementary Schur parameter p^, while 
Gr,r+i and Gt^t+i can contribute each of their entries. In fact, we state the 
following specification of Step E5. 

E5'. Under the assumptions of Step E5, and using the above notations, the 
{i,j) entry of the matrix S equals 



yt, 



where Xr G {ur, Pr, —c^r} and yt G {at, pt, —at} are the entries of 
Gr,r+i and Gt,t+i which can be found as described in the paragraphs 
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abov^ it suffices eacli time to imagine that we are in tlie line segment 
corresponding to tlie current Givens transformation, and tlien imagine 
moving leftwards or rightwards on the path, to obtain the row and the 



column index in ( 30 ) , respectively. 



3.2 Proof of the graphical rule 

The proof that the recipe in Steps E1-E6, E5' leads to the correct form of 
the (i, j) entry of S follows by just expanding the matrix S in an appropriate 
way. Let us sketch here the main steps of the proof. 

PROOF. Throughout the proof, the Givens transformations under the Y\- 
symbol are understood to be multiplied in the order described in Section [l.3[ 
We will assume for definiteness that either i < j or i = j and r < t. Consider 
the given snake-shaped matrix factorization S = YYh=QGk,k+i- Define the 
'sub-snake' 

j 

Sij := Gfc^fc+i- (31) 

k=i-l 

It is clear that the entry of S depends only on the sub-snake Si^j. 
This follows since the other Givens transformations can be considered as 
operations on rows and columns {1, 2, . . . , « — 1} U -|- 1, j -|- 2, . . .} of Stj] 
hence indeed they cannot influence the (?, j) entry of Sij. 

Assume now that si = 1 for some /£{«,... , j}. This means that the line 



segment is positioned to the left of G;_i^;. We can then factor (31 ) as 



W Gk,k+i 1 ' I n Gk,k+i j • (32) 

\k=l J \k=i-l / 



We distinguish between three cases: 



Suppose that / G {i + 1, . . . , j — 1}. The leftmost factor in (32 ) can be 



considered as a row operation acting on rows I, . . . ,j + 1 of the right- 



most factor in (32). By assumption, these row indices are all strictly 



larger than i; hence this factor cannot influence the entry of 



Similarly, the rightmost factor in (32) acts on columns i — 



which by assumption are all strictly smaller than j. We conclude that 



the entry can be influenced by none of the factors G^^k+i in (32), 
and hence it simply equals the entry of the identity matrix, i.e., 
it equals zero. This proves the conclusion in Step E6. 

Suppose that I = i. In contrast to the previous case, we can now only 



conclude that the rightmost factor Gi-i^i in (32) can be removed from 



^Explicitly, Xr is the (i — r, &)th entry of G,.,r+i and yt is the (1 — b,j — t)th entry of 
Gt,t+i, where the boolean b is defined by 6 = if r > t or 6 = 1 if r < t. 
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further consideration. This corresponds to the fact that r equals i 
(and not i — 1) in this case. 



Suppose that I = j. Similarly as in the previous case, we can then 



conclude that the leftmost factor Gjj+i in (32) can be removed from 
further consideration. This corresponds to the fact that t equals j — 1 
(and not j) in this case. 

Getting rid of all the redundant factors G^^k+i as described above, we 
are left with either the identity matrix or with a sequence of Givens transfor- 
mations following a unitary Hessenberg shape (cf . Figure [T| . The relevant 
entries of this matrix can be computed using a straightforward calculation 
and are easily seen to correspond to the given rules in Steps E5 and E5' 
(compare with We omit further details. □ 



3.3 Some corollaries 

A first corollary is the following. 

Corollary 6 (Upper and lower bandwidth of S) The upper bandwidth of 
the snake-shaped matrix factorization S equals the length of the longest sub- 
snake of S whose line segments are linearly aligned in the top left-bottom 
right order (cf. Figure^. Similarly, the lower bandwidth of S equals the 
length of the longest sub-snake of S whose line segments are linearly aligned 
in the top right-bottom left order. 

It follows from Corollary |6] that the unitary five-diagonal matrices C 
have the smallest bandwidth of all snake-shaped matrix factorizations S; 
they have in fact bandwidth 2 in both their lower and upper triangular part 
and hence are five-diagonal. 

A related result on the minimality of the matrix C is the fact [8J that 
any infinite unitary matrix A having lower bandwidth 1 and finite upper 
bandwidth n is 'trivial' in the sense that A can be decomposed as a direct 
sum of matrices of size at most n -\- 1. This result can be shown using only 
some basic linear algebra by noting that under the above conditions on the 
matrix A, this matrix is unitary Hessenberg and hence allows a factorization 
of the form (|4|. The condition on the upper bandwidth of A then easily 
implies that from each tuple of n -|- 1 subsequent Givens transformations 
Gk,k+i in (|4|, there must be at least one for which Gk,k+i has vanishing 
off-diagonal elements; we omit further details. 



A second corollary of the above results can be easily proven from ( 13 ) and 
Lemma [4] Here, the elements of the matrix S given by ( 13 ) are expressed 
in terms of the inner product (11) and the orthonormal Szego polynomials 
(see also Theorem 4.1 in [6j). 
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Corollary 7 By introducing the notation 



Lpi{z) if Si = 0, 
V'iC^) s^ = l, 



then the entries of the snake-shaped matrix factorization S = {111^)1 j>o '^^^ 
given for alii >0 and k > 1 by rji^i = {fi, zfi) and by 



Vi+k,! 



ifi+k, z^+^^+f^fi) if Si+i = ■■■ = Si+k-i = 1, 

other case, 

{fi, fi+k) if Si+i = ■■■ = Si+k-i = 0, 

other case, 



where when k = 1, the condition Sj+i = • • • = Sj+fc-i £ {0, 1} is understood 
to be always valid. 

□ 

As a consequence of Corollary [7] and the graphical rule, by choosing 
appropriate generating sequences one can easily deduce a direct proof of 
Propositions 1.5.8, 1.5.9 and 1.5.10 in |27|. 

4 Connection with Szego quadrature formulas 

In this section we describe some connections between snake-shaped matrix 
factorizations and Szego quadrature formulas. The results in this section 
are actually known for the unitary Hessenberg and unitary five-diagonal 
cases, and the extension to a general snake-shaped matrix factorization S 
turns out to be rather trivial. Nevertheless, we include these results here for 
completeness of the paper. 

Throughout this section, we shall be dealing with a fixed measure fj, as 



described in Section 2.1 and we will be concerned with the computation of 
integrals of the form 

W) := I f{z)di,{z) = r f{e'')dfi{e), 

JT J -TV 

by means of so-called Szego quadrature formulas. Such rules appear as the 
analogue on the unit circle of the Gaussian formulas when dealing with 
estimations of integrals supported over intervals on the real line M. For a 
fixed positive integer n £ N\{0}, an n-point Szego quadrature is of the form 

n 

In{f) ■=^>^jf{zj), Zj^T, j = l,...,n, Zj^Zkiijy^k, 
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where the nodes {zj}'^^^ and weights {Aj}"^^ are determined in such a way 
that the quadrature formulas are exact in subspaces of Laurent polynomials 
whose dimension is as high as possible. The characterizing property is that 
In{L) = I^{L) for all L E spanjz-^ : j = — n + 1, . . . ,n — 1} (the optimal 
subspace): see e.g. [I2l Hg , [201 Chapter 4]. 

In what follows, we will use the notations 7i, C and 5 for the unitary 
Hessenberg, unitary five-diagonal and unitary snake-shaped matrix factor- 
ization induced by the generating sequence {pnjn, respectively. As we have 
already seen, these matrices can all be factorized as HfcLo where the 

Gk,k+i are canonically fixed by ([s]) and ([5]), but where the factors under the 



]^-symbol may occur in a certain order (cf. Section 1.3). 

We start with the following result, which seems to be essentialljUi] due 
to Gragg. It is the unitary analogue of a well-known result for the Jacobi 
matrix when the measure /i is supported on the real line M. 

Theorem 8 ( Gragg 1 18^ ) The eigenvalues of the principal n x n submatrix 
of the unitary Hessenberg matrix 7i are the zeros of the nth Szego polynomial 

Here with the principal n x n submatrix of TC we mean the submatrix 
formed by rows and columns up to n — 1 of 

Proposition 9 (Watkins I33f . Cantero, Moral and Velazquez ^) Theorem 
also holds for the unitary five-diagonal matrix C, i.e., the eigenvalues of 
the principal nxn submatrix of C are the zeros of the nth Szego polynomial 

The above results hold in fact for any snake-shaped matrix factorization 
S (see further). 

For the present discussion, a drawback of Theorem [8] and Proposition |9] 
is that the principal nxn submatrix of {TC, C,S} is in general not unitary 
anymore and hence has eigenvalues strictly inside the unit disk. This means 
that these eigenvalues are not suited as nodes for the construction of an 
n-point Szego quadrature formula. 

The solution to the above drawback is to slightly modify the principal 
nxn submatrix of S in such a way that it becomes unitary. Its eigenvalues 
will then be distinct, exactly on the unit circle T and turn out to be precisely 
the required set of nodes. 

To achieve this in practice, Gragg [T8] and also Watkins [33] introduced 
the idea to redefine the (n — l)th Givens transformation Gn-i,n by 



Gn~l, 



e'" 
e^^" 



(33) 



" Theorem [s] is not explicitly stated in [T^, but it can be easily deduced from the results 
in that paper. 
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where 0, G M denote arbitrary parameters (the second of them will actually 
be irrelevant for what follows). 

With this new choice of Gn-i,n, we can 'absorb' the factors e*^, e*^ in the 
previous and next Givens transformation Gn-2,n-i and Gn,n+i, respectively. 
This means that we redefine 



Gn- 



2,n-l •- 



Gn- 



2,n-l • 



1 

e*^ 



if Sn- 



0, 



(34) 



while in case = 1 we redefine Gn-2,n-i by the same formula (34) but 
now with the factors multiplied in the reverse order. Similarly, we redefine 



Gn.n 



+1 



e"' 
1 



G 



n,n+l ) 



if Sr 



(35) 



while in case — 1 we redefine G^^n-]-! by the same formula (|35[) but now 



with the factors multiplied in the reverse order. We can then put 



G 



n— 1,71 



I. 



(36) 



Note that after the above updates, the value of the snake-shaped matrix 
factorization S remains unchanged but we have succeeded to transform the 
Givens transformation Gn-i,n in (33) into the identity matrix I. Then it 



is easily seen that the snake shape of S can be 'broken' into two pieces, in 
the sense that S = UV where U = llk=o ^k., k+i is the submatrix formed by 
rows and columns 0, . . . , n — 1 of 5, while V = Hfcln Gk,k+i is the submatrix 
formed by rows and columns n, . . . , oo. Note that the matrices U and V have 
a complementary zero pattern and hence they commute with each other. 

Let us now denote with := 11^=0 ^k,k+i the topmost part of the 

'broken' snake S. Note that 5^-1 is a snake-shaped matrix factorization of 
size 71 X n; in particular it is still unitary. Note also that this matrix depends 
on the parameter S M by means of ( 34 ) . 



Remark 10 The principal n x n submatrix of S can be obtained in the 
same way as above, but now replacing the role of e*^ € T in (33) by the 



original matrix entry a^. Note however that lies strictly inside the unit 
disk and hence the resulting nxn submatrix is not unitary anymore; cf. the 
motivation earlier in this section. 



One has then the following result. 

Theorem 11 (Gragg ^18^) Let 9 G R be fixed. Using the above construction, 
the eigenvalues ofTCn-i are distinct, belong to T and appear as nodes in an 
n-point Szego quadrature formula for the measure /i. The corresponding 
quadrature weights are the first components of the normalized eigenvectors 
ofHn-i- 
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Proposition 12 (Watkins 133]) Theorem 11 also holds for the matrix Cn 



Here with 'normalized' eigenvectors we mean that the eigenvectors should 
be scaled in such a way that they form an orthonormal system and that their 
first components are real positive numbers. 

The characteristic polynomial of the above matrix 7in-i (or equivalently, 
Cn-i) is known as a monic para- orthogonal polynomial of degree n ^23^. Note 
that this polynomial depends on the free parameter 9, and hence there is 
in fact a one-parameter family of para-orthogonal polynomials (and so, a 
one-parameter family of n-point Szego quadrature formulas for /i). 

Now one could ask why there is such a similarity between Ti. and C in the 
above results. This is explained by the following basic observation, which is 
essentially due to Ammar, Gragg and Reichel [1] for the case of TC and C. 

Proposition 13 (Based on Ammar, Gragg and Reichel 11]) Let 9 £ M be 
fixed. Then the eigenvalues and the first components of the normalized eigen- 
vectors of Sn~i depend on the Schur parameters hut not on the shape of the 
snake. 

PROOF. Recall that the snake-shaped matrix factorization is given by 
Sn-i = nfc=i ^^^,^+1) for some order of the factors. But it is a general fact 
that the matrices AB and BA have the same eigenvalues; this follows from 
the similarity transformation 

AB ^ A-^{AB)A = BA. (37) 

By applying this idea recursively for the choice A = 11^=^ ^fc,fc+ii for I = 
n — 2, . . . , 1 (only those indices / for which si = 1 have to be treated), one can 
succeed to rearrange the Givens transformations of Sn-i into the unitary 
Hessenberg form Gq^iGi^2 - ■ ■ Gn-2,n-i (compare with ([6|). It follows that 
the eigenvalues of Sn-i are indeed independent of the order of the factors 
Gfc^fc_i_i, i.e., they are independent of the shape of the snake. 

The same argument also shows that the first components of the normal- 
ized eigenvectors are independent of the shape of the snake. To see this, 
consider the eigen-decomposition Sn-i = UDU* , where D is a diagonal ma- 
trix containing the eigenvalues, and f/ is a unitary matrix whose columns 
are the eigenvectors, scaled in such a way that the first row of U has real 
positive entries. The point is now that the only Givens transformation of 
Sn~i acting on the 0th index is Go,i; but in the above argument the latter 



can only appear as the i?-factor in (37), and hence the first row of U is easily 



seen to remain unchanged under the similarity (37). □ 



Corollary 14 Theorems\8\and 11 hold with Ti. replaced by any snake-shaped 
matrix factorization S. 
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PROOF. This follows from Proposition 13 and Remark 10, □ 



Note that Proposition [13] implies that the eigenvalue problems for the 
matrices TLn-i-, Cn-i and Sn-i are conceptually equivalent. Interestingly, 
these problems turn out to be also numerically equivalent since, for reasons 
of efficiency and numerical stability, the eigenvalue computation for {7in-i, 
Cn-i, Sn-i} should preferably be performed using their factorization as a 
product of Givens transformations, rather than using their entry- wise ex- 
pansions. 

Finally, we mention that the development of extensions of Szego quadra- 
ture formulas and the investigation of the connection between them and 
Gauss quadrature formulas on the interval [—1,1] are active areas of re- 
search: see e.g. [3l H H [T2 l [22] and references therein found. A whole variety 
of practical eigenvalue computation algorithms for unitary Hessenberg and 
five-diagonal matrices has already been developed in the literature. In [26], 
Rutishauser designed an LR-iteration. Implicit QR-algorithms for unitary 
Hessenberg matrices were described and analyzed in [HI [lU [191 EH] • In [2 [2T] 
and the references therein, divide and conquer algorithms were constructed. 
Other approaches are an algorithm using two half-size singular value decom- 
positions m , a method involving matrix pencils , and a unitary equivalent 
of the Sturm sequence method [5]. 
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